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Abstract 

The heat trace asymptotics are discussed for operators of Laplace type 
with Dirichlet, Robin, spectral, D/N, and transmittal boundary condi- 
tions. The heat content asymptotics are discussed for operators with time 
dependent coefficients and Dirichlet or Robin boundary conditions. 



1 Introduction 

Standard Dirichlet and Neumann boundary conditions appear in numerous 
physical applications, some of which are nicely described at this Conference. 
In certain cases physics requires consideration of more 'exotic' boundary value 
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problems. For example, divergences of the Casimir energy in non-static, but 
reasonably slow varying, external fields are related to the asymptotics of the 
Schrodinger equation with a time dependent Hamiltonian. After the Wick ro- 
tation the latter are defined by the heat trace asymptotics for operators with 
time dependent coefficients. It is easy to imagine a physical experiment when 
temperature of a part of the surface of a body is kept constant while the heat 
flow from the other part to the outside is negligible. Such physical experiments 
are described by the D/N boundary value problem. Transmittal boundary con- 
ditions appear in the case of semi-transparent surfaces or when the geometry 
of the manifold is not smooth. The most fashionable example (and the closest 
to the topic of the present Conference) of the non-smooth geometries is given 
by the brane world scenario p3| . Spectral boundary conditions are of relevance 
in one- loop quantum cosmology and supergravity jl5], Furthermore, given 
their nice transformation properties under chiral rotations and supertranslations 
there is little doubt that study of spectral boundary conditions is also useful. 

Let D be an operator of Laplace type acting on the space of smooth sections 
to a vector bundle V over a compact Riemannian manifold M of dimension m 
with smooth boundary dM. Let Dg be the realization of D with Dirichlet or 
Robin boundary conditions; we will consider more 'exotic' boundary conditions 
presently. Let e~ tDs be the fundamental solution of the heat equation; u := 
e~ tDe 4> is determined by the equations: 

u{x; 0) = cf>, Bu = 0, and (d t + D)u = 0. 

Let / be a smooth localizing or smearing function. We define the smeared 
heat trace: 

a(f;D,B)(t) := Tr L2 (/ e ~ tDB ). 
As t I there is an asymptotic series |^5|, Q 

a(f; D, B) ~ £„>o ^"""^M/, A B). (1) 

The asymptotic heat trace coefficients may be decomposed as the sum of an 
interior and a boundary contribution: 

a n (.f, D, B) = af (/, D) + aP n M {f, D, B). 

The invariants a^f and a® M are computable as integrals of local geometric in- 
variants. 

Let ip be an auxiliary section to V defined over the boundary dM and let the 
potential p measure internal heat sources and sinks. Let u be the temperature 
distribution defined by the inhomogeneous equations: 

u(x; 0) = 4>, Bu = ip, and (dt + D)u = p. 

With Dirichlet boundary conditions, we keep the boundary at constant tem- 
perature ip' with Neumann boundary conditions, we pump heat into M at a 
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rate defined by ip to control the heat flow in the normal direction. Let p be the 
specific heat of the manifold. The total heat content 



P{pA^,P\D,B)(t) := J M up 
has an asymptotic expansion as t { 

The heat content asymptotics (3 n can be decomposed as the sum of an interior 
and a boundary contribution given by locally computable invariants. 

The coefficients a n and p n encode spectral information about the global 
geometry of the manifold. In Section H we discuss the interior invariants azf. 




These invariants vanish if n is odd. In Theorem 2.1, we give formulas 
the invariants a n for n = 0, 2, 4; formulas for the invariants ag fisfl and 
are known. 

In Section pi we define the Dirichlet and Robin boundary operator - see 
equation (^). In Theorem [O, we give formulas |9|, ^?], [Io| for the associated 
boundary correction terms a^ 1 if n < 4; formulas for 05 are known Oj. 

In Sectio n p[ we define transmittal boundary conditions - see equation (^|). 
In Theorem |4.l| , we give formulas ^ p3|, p9| for the boundary correction terms 
O'n if n < 3; formulas for a a are knownp3[7 

In Section ^, we discuss spectral boundary conditions. In contrast to Dirich- 
let, Robin and transmittal boundary conditions, spectral boundary conditions 



are non-local. In Theorem 5.1, we give formulas for the boundary correction 
terms if n < 3 Jl3| , pl[ . Apart from normalizing factors involving powers of 



47T, the formulas of Theorems 2.1, 3.1, and 4.1 involve coefficients which are 



independent of the di men sion m of the underlying manifold. In contrast, the 



formulae of Theorem 5J are very dimension dependent. This is one of the 
notable features of spectral boundary conditions. 

In Section ||, we consider a time dependent family V = T> t of operators of 
Laplace type. The heat temperature distribution is defined by: 

u(x; 0) = 0, Bu = 0, and (d t + V t )u = 0. 

The map <f> — > u is described by a smooth kernel function K with the property 
that: 

u(x; t) = J M K(t, x, y, V, B)(j)(y)dy. 

The heat trace asymptotics are then defined not by the heat trace but directly 
in terms of the kernel function: 

a(f, V, B)(t) := J M /Tr K(t, x, x, V, B) 
~£ n t (n - m)/2 a„ 
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In Theorem 6.1 we give formulas for the interior invariants. We define boundary 
conditions in equation (Q) which are time dependent. In Theorem 6.2, we give 
formulas for the boundary correction in the heat equation asymptotics. 

In Section |?], we give a brief discussion of the D/N problem [l2|, Here, 
in contrast to other boundary conditions, there is not a classical asymptotic 
expansion at the et 3 level. 



In Section pL we discuss the heat content asymptotics. In Theorem ST 
we give formulae || [?], |24| for the invariants f3 n for n < 4 for Dirichlet or 
Robin boundary conditions. The coefficients which appear do not depend on 
the dimension m. In the static setting, partial results are available for fa flq, p|. 



2 Interior Heat Trace Coefficients 

We introduce the following notational conventions to describe the interior heat 
trace coefficients . Let Greek indices fi, v range from 1 to m and index a 
local coordinate frame. Let Latin indices i,j,k 7 l range from 1 to m and index 
an orthonormal frame. We adopt the Einstein convention and sum over repeated 
indices. The operator D determines a connection V and an endomorphism E 
so that 

D = -(TrV 2 +E). 

If we express 

D = -(g» v d fl d I/ +a»d fl + b) 

relative to a local system of coordinates, then the connection 1 form to and the 
endomorphism E are given by: 

w s = \g v s(a v + gi u 'T lur v I)> and 

E =b- g Vfl (d u w fJ , + uj u LUft - uj a T vl f). 

If D = 8d is the scalar Laplacian, then the connection V is trivial and the 
endomorphism E vanishes. More generally, if D = (dS + 5d) p is the Laplacian on 
p forms, then V is the Levi-Civita connection and E is given by the Weitzenbock 
formulas |fl~s|] . If D is the spin Laplacian, then V is the spin connection and with 
our sign convention E = —\t where r is the scalar curvature. 

Let ';' denote multiple covariant differentiation with respect to the connec- 
tion on V and the Levi-Civita connection of M. Let 

Pij ■= Rikk] and r := p rl 

be the Ricci tensor and the scalar curvature. Let Q be the curvature of the 
connection V. If A is a scalar invariant, we let Tr („4) := Tr (-4/). 

The invariants vanish if n is odd. If n is even and if n < 4, then we have 




Theorem 2.1 
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1. a M (/,C) = (4^)-™/2| M Tr(/). 

2. af (/, D) = (4^)-™/2 i | m JTr (r + 6E) , 

3. af (/, D) = (4^)-™/ 2 ^ f /Tr {60£ ;fefe + 60tE + 180E 2 

• :!()<!,,<),, + l2 T . kk + 5r 2 - 2\p\ 2 + 2\R\ 2 }. 

3 Heat Trace Asymptotics for Robin and Dirich- 
let Boundary Conditions 

Near the boundary, let Roman indices a, b range from 1 to m — 1 and index a 
local orthonormal frame {e a } for the tangent bundle of dM. We let e m be the 
inward unit normal. Let 

L a b : = (V ea e;,,e m ) 

be the second fundamental form. Decompose the boundary of M as the disjoint 
union of two closed (possibly empty) sets: 

dM = C N UC D . 

Let u- m be the covariant derivative of u with respect to the inward unit normal 
using the natural connection defined by D. Let the boundary operator 

Bu := u\ Cd © (u- m + Su)\ Cn (2) 

define Dirichlet boundary conditions on Cd and Robin boundary conditions 
on Cat. Let ':' denote multiple covariant differentiation of tensors defined on 
dM with respect to the connection on V and the Levi-Civita connection of the 
boundary Note that ';' and ':' differ by the second fundamental form. We have 

Theorem 3.1 

1. a$ M (f,D,B)=0. 

2. af M (f, D, B) = -(4^-™)/ 2 i J Cd Tr (/) + {A-k)^-^\ Tr (/). 

3. al M (f,D,B) = (47r)-™/ 2 i f Cn Tr {2fL aa - 3/ ;m } 

+ (47r)-™/ 2 i J Cn Tr {2fL aa + 12 fS + 3/ ;m }. 

I 4 M {f, D, B) = -(47r)( 1 -™)/ 2 ^ J Cd Tr {f(96E + 16r - 8/w + 7L aa L bb 
-10L ab L afc )-30/ ;m L aa +24/ ;mm }+(4^)( 1 - m )/ 2 3| 5 J Cn Tr {/(96£+16t 
-8p mm + 13L aa L bb + 2L ab L ab + 96SX aa + 192S* 2 ) + f m {6L aa + 965) 
+24/ ;mm }. 
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5. a d 4 M (f,D,B) = {4ny m / 2 ^ J Cd Tr {/(-120£ ;m + l20EL aa - 18r ;m 

mm-Lbb ^^Rambm-^ab 

— ^fL a bL a bL cc + ^LabLbcLac) + / ;m (— 180i? — 30r — ^-L aa Lbb 
+ fL ab L ab ) + 24f. mm L aa - 30/ ;l4m } + ^tt)-™/ 2 ^ Tr {/(240£ ;m 
+120£ J L aa + 42r ;m + 24L aa:bb + 20ri aa - Ap mrn L bb 12 Rambm L ab 
+4:RabcbLac+^§-L aa LbbL cc +8L a bL a bL cc +^-L a bLb c L ac +720SE+120ST 
+UASL aa L bb + ^SL ab L ab + 480S 2 L aa + 480S 3 + 120S :aa ) + / ;m (180£ 
+72SL QQ + 240S 2 + 30r + YlL aa L hh + 12L ab L ab ) + / ;mro (1205 + 24L aa ) 



4 Transmittal boundary conditions 

Let dM be empty. We suppose given a hypersurface £ which divides M into 
two smooth components M . We also suppose given operators of Laplace type 
D ± on Af ± . Let ^ be the inward normal of E C M + . For <p = we 
define: 

^:={0+|s-r|s} (3) 

®{(v^+)| s -(v,;r)|E-5^ + |E}. 

Thus satisfies transmittal boundary conditions if 4> is continuous and if the 
normal derivatives of <fi + match to the normal derivatives of 0~ modulo the 
impedance transmission term S. We let D& be the realization of D = (D + , D~ ) 
with these boundary conditions. Let / = (/ + ,/~) be smooth on and 
continuous on S; we impose no matching condition on the normal derivatives 
of /. Let 

a(f,D,B) = Tr L 2(fe- tD °) 

We can decompose the invariants a n in the form: 

a n (f, D, B) = eg + (/, D) + a™~ (/, D) 
+a^(f,D,B). 



The invariants a^f can be computed using the formulas of Theorem 2.1. Let 
v^- be the inward normals of S C M^; v = i> + = —v~. We let 

"a ■= V+ - V- and L± := ±(V± e 6 , i/). 

The tensor cj a measures the failure of the connections to agree on S; it 
is a chiral tensor - if we interchange the roles of + and — , then this tensor 
changes sign. The tensors are the second fundamental forms of the inclusions 
E C . We refer to [^3| for the proof of the following theorem; see also related 
work in [|||. 
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Theorem 4.1 

1. a$(f,D,E)=0. 

2. af(f,D,E)=0. 

3. a$(f, D, S) = (47r)-W2i J s T r {2f(L+ a + L~ a ) - 6/5}. 

* a? (/, D, S) = (4 7 r)( 1 — )/ 2 ^ J E Tr L+ + L" L£ + 2L+ L 6 " 6 ) 

+3/(^+ 6 + + 2L+ + 9(L+ + L- )(/+ + + /"_) 

+48/S 2 + 24/a, u; a - 24/(L+ + L" )S - 24(/+ + + /"_)H}. 

5 Spectral boundary conditions 

Let P : C ca (Ei) — » C X, (E 2 ) be an elliptic complex of Dirac type; this means that 
the associated second order operators P*P and PP* are of Laplace type. Since 
such an elliptic complex does not in general admit local boundary conditions 
H , we impose spectral boundary conditions. Let 7 be the leading symbol of the 
operator P. Then 7 + 7* defines a unitary Clifford module structure on E\®E2- 
Let V = Vi © V2 be a compatible unitary connection ju|. This means that: 

V(7 + 7*) = 0, and 
(Vs,s) + (s, Vs) = d(s,s). 

In general this auxiliary connection will not coincide with the connections asso- 
ciated to the Laplacians Ai = P*P and A2 = PP*. We expand 

p = + v 

where ip is a smooth linear map from Pi to p2- We parallel translate frames for 
E along the normal geodesic rays defined by the inward unit normal. Relative 
to such a gauge, we have V m = d m . We set x rn = to define the tangential 
operator B on C°°{E 1 \qm)- 

B-= (7 m )^{7 Q V Q +^}. 

Let B* be the adjoint of B relative to the structures on the boundary and let 
Q be an auxiliary self-adjoint operator. We define 



Let B denote projection on the span of the eigenspaces corresponding to the 
non-negative spectrum of A. Let Pg be the associated realization of P and 
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let Db '■= {Pb)*Pb- Results of Grubb and Seeley show that there is an 
asymptotic series as t j of the form: 

Tr L2 {fe- tD *} 
~ Zo<k< m -iak(F,D,B)t( k -™y* 

modulo terms which arc 0{t~s). (There is in fact a complete asymptotic series 
with log terms, but we shall only be interested in the first few terms in the 
series). We shall assume m > 4 so the terms a n for n < 3 are well defined. Let 



c(m) :=r(f )r(i)- 1 r(2s±i)- 1 . 



V 1 := 7m V>> and 

2 I'-y-il V 2 

We refer to jl3|, |l| for the proof of the following theorem: 
Theorem 5.1 We have 

1. a (f,D,B) = (^)- m / 2 J M Tr{f}. 

2. ai (f,D,B) = (47r)-(™- 1 )/ 2 i(C(m) - 1) Tr {/}. 

3. a 2 (f, D,B) = af (/, £>) + (4tt)-W2 J SM /Tr + </>*] 

+ |(1 - |^C(m))L aa } - 3^(1 - ±7rC(m))Tr{/ ;m }. 

^. a 3 (f,D,B) = (47r)-(— 1 )/ 2 / aM /Tr{3L(l - 2M)(^ + ^*) 

+ ^ (5 - 2m + 7 - 8 ^+ 2 2m2 C(m) )^* + m ± = ^ (2m - 3 - 2ro2 ~^+ 5 C(m) ) 
•(7lW + 7^*7^*) + ig^ny (1 + ^ C(m)) 7a T ^* 

i L ab L ab / Vr±5m , 23-2m-4m 2 r'fm'l') i L aa L bb ( 17+7m 2 
+ 48(m+l) 1 - 4 "r m-2 Kj V' 1 ')) ~T 48(m 2 -l) \ 8 

+ V-lW + 5,n- l c(m)) + _^_ C(m)(ee + _L_ 7 T0 7 T e)} 

+ fe^t(^ - ^C(m))Tr{/} + I ^l gj (2C(m) - l)/ ;mm Tr {/}. 

6 Time dependent coefficients 

Previously, we have considered static operators. Let T> be an operator of Laplace 
type where the coefficients are time dependent. We expand 

Vu := Du + Y^r>0 tT {Gr,i 3 U;ij + ^r,iU-i + £ r u] 

and consider time dependent Dirichlct and Robin boundary conditions: 

Bu := u\ Cd © (u- m + Su + t(T a u- a + Siu))\c N - (4) 
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We consider the heat equation: 

u(x; 0) =4>, Bu = 0, (d t + V)u = 0. 
There is a smooth kernel function K(t, x, y, 2?, £>) so that we may express: 

u(x;t) = f M K(t,x,y,V,B)(j)(y). 
We take the fiber trace to define 

a(f,V,B) = J M fTi v K(t,x,x,V,B) 

~ £„> * (n - m)/2 M/,^£)- 

This agrees with the previous definition if V is static. We refer to |^| for 
the proof of the following two results which give the additional terms in the 
asymptotic expansion arising from the time dependent nature of the coefficients: 

Theorem 6.1 

1. a™{f,V) = a™{f,D). 

2. af{f, V) = af{f, D) + (4tt)— ' 2 \ f M /Tr 

3. of (/, V) = af (/, £>) + (4^)- 1 / 2 ^ J M /Tr (fS,^,^ + f G Uj G Uj 

+6OC/2. a — 180fi + l5Gi,uRjkkj — 30Gi,ijRikkj 
+90Gi,uE + + ir,c; ; ,,,, - :Wr. ,,:,,)■ 

Let So denote the associated static boundary conditions. We have: 
Theorem 6.2 

1. a d n M (f,V,B) = a^ M {f,D,B ) for n < 2. 

2. a d 3 M (f,V,B) = al M (f,D,B ) + (4tt)( 1 — V 2 ^ /Tr (-240 llOO ) 

+ (47r )(i-™)/^/ Cjv /Tr(240 1 , aa ). 

3. af M (/, £>, B) = af M (/, D, B ) + (4tt)-™/ 2 JL / c d /Tr {3O0 1)aa L 66 

-6001 

.mm ,mm;m 

3O0i 

,aa;m 

+ O0i 

,am;a 

30^i, m } 

+/ ;m Tr {-450i,a« + 450 x , mm } + (4^)- m /2_i_ J Cn /Tr {3O0 Ma L 66 

+ 1 20Gl : mmLbb — 15001, a bi a fc — 6001, mm;m + 6001, aa;m + 150.Fl, m 

+18O50i, aa - 18O50i, mm + 3605i} + / ;m Tr {450i, aa - 450i, mm }. 



7 The D/N Problem 



In Section pL we assumed that Cjv H Co was empty to define the boundary oper- 
ator B of equation (Q) . This meant that the Neumann and Dirichlet components 
did not interact. In this section, we suppose that S :— Cd HCjv is a non-empty 
smooth submanifold of dM of dimension m — 2. 

We can motivate this more generalized setting with a physical example. Let 
M be a solid ball which floats in ice water. The part of the boundary of the ball 
which is in air satisfies Neumann conditions and the part which is in the water 
satisfies Dirichlet conditions. Thus B is defined by complementary spherical 
caps about the north and south poles of the ball which intersect in a circle of 
latitude. 

The setting where £ is not empty is known in the literature as the N/D 
problem. It has been investigated extensively from the functional analytic point 
of view |2^, ||^, ^5|. It is natural to conjecture the asymptotic expansion 
described in (fiT) could be generalized to this setting by adding an extra integral 
over £ of some suitably chosen local invariant. Preliminary computations [l2| 
suggest the additional correction term for n — 2 is given by: 

«2 =-f(4^)- m/2 / E Tr(/). 

However, it has been shown |L4| that the asymptotic expansion does not ex- 
ist with locally computable coefficients at the 03 level. Thus probably either 
log terms arise or non-local terms arise; it is also possible, of course, that no 
asymptotic expansion exists. 



8 Heat Content Asymptotics 

Let I? be a time dependent family of operators of Laplace type. Let ip(y> be a 
smooth section to V defined over dM. On the Neumann boundary component 
C/v, we use a Neumann heat pump to pump heat into M at a rate defined by tp; 
in this setting, the parameter S controls the coupling between the heat transfer 
and the temperature difference on the Neumann component. On the Dirichlet 
component we use a Dirichlet heat pump to keep the temperature at ip. Let p 
be a heat source. The temperature distribution u = u p ^^(x; t) which is defined 
by these data is the solution to the equations: 

(dt + V)u = p, u(x; 0) ~ (/), and 
Bu = tp. 

Let p be the specific heat; we regard p as a section to the dual bundle V* and let 
(•, •) denote the dual pairing between V and V* . The total heat energy content 
P is defined by (3(t) := J M up. We expand in an asyptotic series as t J, to 
define the associated heat content asymptotics: 

P~En tn/2 Pn(p,<f>^,P;D,B). 
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Let D and B be the dual operator and dual boundary condition on the dual 
bundle V*. We summarize results of §, g, |(], 

Theorem 8.1 

1. p (p, </>,,/;, p;V,B) =J M {<t>,p). 

2. 0i(p, 0, V, p; V, B) = -^J CD {(<f> - Vo, p)}. 

3. p 2 (p,<l),iP,p;V,B) = -J M {(D4>,p) - (p ,p)} 

+/ Cjv {<(B0-<M,P>}. 
* fofah^fr^B) = -^Ic D ii(Po,p) - I(D^p) 

-j{(<f>-lfa),pp) + 5<(0-^):a,P:o) ~ |(^1,P> + {{~\E + ±L aa L bb 
— gL a bL a b + -gRamam ~ Gl,mm) ' {4> ~ "0o), P)} 

+ s^!c N {W-^),Bp)}- 

5. fa(p,<i>^,p;V,B) = \j M {(pup) - (D Po ,p) + (D^Dp) 

-((Gl,ij<f>;ij +Fl,i<t>;i +£l4>),p)} + S c D il L aa(po, p) ~ l(P0,P;m) 
-zLaa(lh,p) + WuP;m) + P) + J<(0 ~ M (£p);m> 

-^{LaaDcf), p) - i(L aa (0 - ip ),Dp) + {{\E. m - ±L ab L ab L cc 

~^"^L a i ) L ac L/i )C — YQ-^ambm^ab ~l~ Jg Rabcb L ac ~f~ S2^~' Tn 
+ j^£ a b:ab)0 ~ Vo), P) ~ \ L ab{Q> ~ 1po):a, P:b) 
-g(^amO ~ 1p0):a,p) + | (^am (</> ~ ^o),P:a) 
+ (^5l,mm;ra ~ jGl,mmL a a ~ Tg-^l.ro) ((^ ~ Ipo) ■> P) 
am(0 - V'oXa, P) + ^Gl, mm {{<P ~ Ipo), P;m)} 

+J CN {l(Bp ,p) - \{{Bcf>- fl> ),Dp) - l(D<p,Bp) - \(^ P ) 
+{{\S + \L aa ){Bcl> - fo), Bp) - \Gi, mm {{B<t> ~ V>o), p)}. 

9 Remarks 

We have presented explicit combinatorial formulas for both the heat content and 
the heat trace asymptotics. One of our motivations in computing these invari- 
ants was to see if there was a direct combinatorial link between the invariants; 
there does not seem to be one immediately evident although techniques used 
in the computation of both the heat content and the heat trace asymptotics 
share certain common features and in principle there are methods which might 
permit both to be computed simultaneously. Another example of an asymptotic 
formula involving geometric data arises from expanding the volume of a tube 
of radius r around a submanifold N embedded in an ambient manifold, see for 
example |3(| . Again, there does not seem to be any direct combinatorial link 
between these asymptotic formulae and those we have presented here. 
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